Tunable Electron Interactions and Fractional Quantum Hall States in Graphene 
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The recent discovery of fractional quantum Hall states in graphene raises the question of whether 
the physics of graphene offers any advantages over GaAs-based materials in exploring strongly- 
correlated states of two-dimensional electrons. Here we propose a method to continuously tune 
the effective electron interactions in graphene and its bilayer by the dielectric environment of the 
sample. Using this method, the charge gaps of prominent FQH states, including v = 1/3 or u = 5/2 
states, can be increased several times, or reduced to zero. The tunability of the interactions can be 
used to realize and stabilize various strongly correlated phases and explore the transitions between 
them. 

PACS numbers: 63.22.-m, 87.10.-c,63.20.Pw 



Introduction. Two-dimensional electron systems 
(2DES) placed in a high magnetic field exhibit strongly 
correlated phases characterized by fractionally quantized 
Hall conductivity ll| , quasiparticles that carry a fraction 
of electron charge [3 , and fractional ( Abelian or possibly 
non-Abelian) statistics [H, Q . These remarkable phenom- 
ena occur in the extreme quantum limit - the fractional 
quantum Hall (FQH) regime - when the number of elec- 
trons, TVe, is comparable to the number of magnetic flux 
quanta through the 2DES, N^, corresponding to a partial 
filling ly = Ne/N^ of one of the lower Landau levels (LLs). 
When v is swept through the series of simple fractions in 
the lowest ti = Landau level (LLL) where Ne < N^, the 
electrons condense into Laughlin states which describe 
the strongest observed fractions (i^ = 1/3, 1/5, . . .) m, or 
weaker states that belong to the so-called hierarchy [J, 
or composite fermion series Within the n ~ 1 LL, 
the effective interaction changes due to the nodal struc- 
ture of LL orbitals, and some of the more exotic states, 
such as the Read-Rezayi states 0], are more likely to be 
favored. The experimentally most important member of 
the Read-Rezayi series is the Moore-Read (MR) "Pfaf- 
fian" state Isl, believed to describe the FQH plateau at 
v = 5/2 [1, HI . Quasiparticles of the Read-Rezayi series 
obey the non-Abelian statistics 0] which is of interest for 
topological quantum computation (loj . 

Although many FQH states have been discovered in 
GaAs-based 2DES, these systems are plagued by the fact 
that their 2DES is buried inside a larger 3D structure. 
This fixes the effective interactions at values that arc of- 
ten not optimal for some of the most interesting FQH 
states, including the Read-Rezayi series. Theoretically, 
such states are known to be verysensitive to the form 
of the effective interactions [111, ■ Another problem 
stems from the strong dielectric screening and finite well- 
width [l3| in GaAs, which weaken the electron-electron 
interactions, thereby making FQH states fragile. This 
has been a major obstacle in the studies of the possi- 
bly non-Abelian states, which could only be observed in 



ultra-high- mobility samples [§] (see however Thus, 
it is desirable to find an alternative high-mobility 2DES 
with strong effective Coulomb interactions that are ad- 
justable in a broad range. 

A promising candidate for this kind of material is 
monolayer graphene (MG), a high-mobility atomi cally 
thick 2DES [15[, where recently a v = 1/3 state [iq . 
and several other states [l^l, have been discovered. Re- 
markably, due to graphcne's truly 2D nature, the short- 
range electron interactions greatly exceed those in GaAs, 
which leads to a significantly more robust FQH state 
at J/ = 1/3 Il8|. A closely related material, bilayer 
graphene (EG) [l5| has similarly high mobility, and ex- 
hibits interaction-induced quantum Hall states at integer 
filling factors at low magnetic fields [31 . This indicates 
that, similarly to MG, the underlying electron interac- 
tions are strong and one could expect robust FQH states 
in EG as well. 

Main results. Here we propose a way to continuously 
tunc the interactions in graphene in a wide range, using 
the fact that 2DES in graphene is exposed and its prop- 
erties can be controlled by the dielectric environment, 
as illustrated in Fig. [TJ Using exact diagonalization cal- 
culations, we show that the tunability can be used to 
significantly increase the excitation gap of = 1/3 state 
in both MG and BG, as well as in the half-filled n — 1 LL 
in BG. The latter, similarly to the case of GaAs, resides 
in the gapped MR PfafRan state 11 , • As the gap 



is varied, the overlaps between the exact states and the 
model wavefunctions improve by a few percent, and their 
topological character becomes better protected. The re- 
duction of the gap induces the transitions to compressible 
and crystalline phases |12| . 

Model. A graphene sample is situated in a dielectric 
medium with permittivity ei, and a semi-infinite dielec- 
tric plate with permittivity £2 7^ ci is placed at a dis- 
tance d/2 away from the graphene sheet. The effective 
interactions between electrons in graphene change due 
to the surface charges induced at the boundary between 
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FIG. 1: (Color online). Setup of a grapliene system with 
tunable electron interactions. Electron states are affected by 
a dielectric plate placed in the vicinity of the surface. By 
varying the dielectric permittivity of the plate and its distance 
d/2 to the graphere layer, interactions of the desired form can 
be engineered. 
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Below we measure the distance d in units of the magnetic 
length Ib and the energy in units of e^/eils- Ratio d/ls 
controls the effective interactions within a partially filled 
LL (see below). All the gaps quoted here should be mul- 
tiplied by a factor ecaAs/ei if comparison is to be made 
with GaAs 2DES. An important advantage of this setup 
is that the interactions can be tuned in situ by varying 
the magnetic field B, which modifies the ratio d/ls- 

We project the interaction to a partially filled LL under 
consideration, following the standard procedure in FQH 
numerical studies [1^. Within a LL, V{r) is parameter- 
ized by the Haldane pseudopotentials {Vm\ [20[, which 
can be conveniently evaluated from the Fourier transform 
of the Coulomb interaction, V^(q), and the form factor 
F(q) encoding the properties of LL orbitals. Due to the 
difference in band structure, the form factors in MG, BG 
are generally distinct from those of GaAs [2lj . 

For what follows, it is useful to review the previous 
results on the relation between pseudopotential values 
and the stability oi v = 1/3 and v = 5/2 states. At 
V = 1/3, the bare Coulomb interaction in n = LL fa- 
vors the Laughlin state. Reducing Vi while keeping Kra>3 
constant eventually destroys the gap and a compressible 
state sets At half filling of the LLL {v =1/2), the 

Coulomb interaction gives rise to a gapless Fermi liquid 



of composite fermions [22] ■ However, in half- filled n = 1 



LL there is a fully developed plateau in experiments [8| , 
attributed to the MR state [1] . MR state is an eigenstate 
of a particular three-body repulsive interaction [9| . Re- 
markably, in numerical studies, the ground state of the 
Coulomb interaction at = 5/2 is seen to be adiabat- 



ically connected to the MR state [23|, [2^; the overlap 



of the ground state with the MR state is improved by 



FIG. 2: (Color online). Ratio V1/V3 as a function of a for 
several values of d/ls in n = LL of MG and BG. Inset: 
typical ratio V1/V3 for a GaAs infinite quantum well of width 
w/Ib, plotted on the same scale. 



the increase in Vi pseudopotential (or, alternatively, by 
reducing V3) [13] ■ Therefore, theory shows that varying 
the first few pseudopotentials provides a convenient way 
to assess stability and induce transitions between FQH 
states. However, so far it has been difficult to find a con- 
trolled way of tuning Vm 's experimentally in a sufficiently 
broad range. 

In Fig. [2]wc show an example of the ratio V1/V3 cal- 
culated for n = LL of MG and BG using the setup 
illustrated in Fig. [TJ Changing a and d/ls leads to a 
large variation of V1/V3 with respect to its pristine value 
(a = 0). For comparison, wc also show a typical vari- 
ation of V1/V3 achievable in GaAs wide quantum wells 
(inset). In what follows, we investigate the model de- 
fined by Eq. ^ using exact diagonalization studies [2^ . 
We focus on the case of spin- and valley-polarized FQH 
states, such as the Laughlin or MR state, and consider 
the physics of a single non-degenerate Landau sub-level. 
For example, in case of the MR state, wc expect the to- 
tal filling factor i/ = -4 + 2M + 3/2 of the BG to be the 
most suitable for realizing the Pfaffian state. Note that 
the exchange interactions favor the following order of fill- 
ing the eight n = 0, 1 Landau sub- levels of the BG [25| : 
|0(TiSi), |l(TiSi), \Ocr2S2), |lcr2S2), where lats^), 1 < z < 4 
are the four orthogonal states in the space of internal 
indices, with s and cr denoting two spin and two valley 
species, respectively. At ly = —4 + 271/ + 3/2, where 
M = 0,1,2,3, the M pairs of 71 = 0,1 sub-levels are 
expected to be filled, and in the top-most pair, n = 
Landau sub-level is filled completely, as dictated by the 
exchange interactions, while n ~ I sub-level is half-filled. 
For brevity, below we refer to this n = 1 sub-level as 
"n = 1 LL in BG" . 

One important criterion of the stability of FQH states 
is their gap to excitations, which can be either neutral 
(quasiparticle-quasihole pair above the ground state) or 
charged (created by changing the magnetic flux by one 
unit). These two gaps are generally different, but for 
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FIG. 3: (Color online). Charge gaps for A'^ = 10 particle 
V = 1/3 state in the n = LL of MG and BG (a), and 
iV = 14 particle v = Ijl filled n = 1 LL of BG (b). By 
changing a and djls, the gaps can be increased several times 
or brought down to zero, following similar trends in each of 
the LLs. 



the model defined by Eq. ([T]) we find that they show 
the same behavior as a function of parameters a and 
djlB- We therefore focus on the charge gap, defined by 
Ac = Eqh + Eqy — 2Eo, where Eq is the ground state 
energy and Eqh , Egp are the energies at iV$ + 1 , — 1 , 
respectively {N^ is kept fixed). In evaluating Ac, it is 
convenient to place the FQH system on the surface of 
a sphere, where incompressible ground states are easily 
identified by their zero angular momentum. On a finite 
sphere, FQH states are further characterized by the shift 
S [2^, which relates iVg, Nq, to the filling factor ly in the 
thermodynamic limit via Ng, = N^/i' + S. In Fig. |3] we 
plot Ac for the case of v = 1/3 filling in the n = LL 
(MG, BG) and v=l/2 filling of n = 1 LL in BG. We do 
not show the results for the n = 1 LL of MG because they 
support significantly weaker pairing correlations than in 
BG [23|. In these plots, the system size is fixed at iV = 10 
{ly = 1/3) and N = 14 {i/ ^ 1/2), and S is chosen to be 
—3. This shift corresponds to the case of bare Coulomb 
interaction where the states are known to be described 
by Laughlin and MR wavefunctions. We estimate that 
the gaps can be increased 2 — 3 times with respect to 
the vacuum value (a = 0). The gaps in n = 1 LL of 
BG are overall significantly smaller than those in n = 
LL, which is consistent with the experiments detecting 
fewer FQH states in n = 1 LL of GaAs. When a < 0, 
the gaps become very small and our analysis would likely 
need to be extended to other shifts that may be preferred 
by these compressible states in minimizing their ground- 
state energy. This analysis will be provided elsewhere. 
We emphasize that all the systems smaller than those in 
Fig. [3] show qualitatively identical plots with the same 
relative variation of the gap (we use the rescaled mag- 
netic length and background charge correction in Fig. [3] 
in order to minimize the finite-size effects [fl]). 

Comparing the gaps in Fig. [3] with the ratio V1/V3, we 
note that the maximum value of V1/V3 does not coincide 
with the maximum of the gap. This is because the gaps 
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FIG. 4: (Color online). Overlaps between the exact ground 
state, for each a and d/ls, and the MR Pfaffian state, for 
A'^ = 16 particles on the sphere (a), and = 14 particles 
on the torus (b). Maximum overlap occurs for intermediate 

d/lg. 

are determined by the complete set of pseudopotentials - 
while negative a's increase the ratio V1/V3, they also re- 
duce the magnitude of each individual Vm , thereby lower- 
ing the gap. However, the Laughlin states and higher or- 
der hierarchy states are excellent trial states for positive, 
as well as negative, a i.e. for long-range Coulomb-like in- 
teraction as well as strong short-ranged repulsion. This 
is reflected in the high value of the overlap between an 
exact state and the Laughlin wavefunction which varies 
by as little as ±1% over the entire phase diagram in Fig. [3] 
(not shown). 

The stability of the MR Pfaffian state is more subtle 
and we expect sharper variation of the overlap as a func- 
tion of a and d/ls- In Fig.|3]we show the overlap between 
the exact state of the dielectric model (for given a and 
d/ls) and the MR Pfaffian state. The overlaps are re- 
markably consistent for the two choices of boundary con- 
ditions: on the sphere with S = —3, Fig.lD^a), and on the 
torus with the hexagonal unit cell. Fig. |3Jb). Although 
MR state is threefold degenerate and characterize d by 
different pseudomomenta k^; i = 1, 2, 3 on the torus [12l |. 
the sixfold symmetry of the Bravais lattice guarantees 
the degeneracy of these subspaces and makes it possi- 
ble to calculate the overlap in the usual way. The max- 
imum of the overlap occurs for d between Ig and 21b, 
and is somewhat higher on the sphere, due to the bias of 
the shift which lifts the particle-hole symmetry (the anti- 
Pfaffian [i^l, having 5 = 1, occurs in a different Hilbert 
space and does not "mix" with the Pfaffian, like on the 
torus). 

We now present additional evidence in favor of the 
state being the MR Pfaffian (Fig. [5]). First, we track the 
evolution of the state upon further variation of Vi pseu- 
dopotential. Fig. [5l^a). The increase of Vi would partly 
mimic the effects of LL mixing [2^ that we neglected 
so far (we do not investigate the effect of three-body 
interaction that also arises in this case). In Fig. ^a.) 
we show the overlaps between the exact state, |Diel), at 
a = l,d ^ 1.51b (chosen from the high-overlap region in 
Fig-HUb)), the bare Coulomb n = 1 LL state, |C, n = 1), 
with the MR Pfaffian |MR) and its particle-hole sym- 
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FIG. 5: (Color online). (a) Overlaps between the 
(un)screened Coulomb interaction eigenstate and MR Pfaf- 
fian, as well as particle-hole symmetrized MR, for = 14 
particles on the torus. Screened Coulomb eigenstate is de- 
fined by a = l,d = 1.51b (compare with Fig. Iljb)). (b) 
Entanglement spectrum on the sphere for A'^ = 16 electrons, 
for bare and screened Coulomb interaction. Entanglement cut 
is chosen as 0|1 in the notation of Ref. [s^; other cuts yield 
similar results. 



permittivities and thicknesses of the order of Ib placed in 
the vicinity of the surface. In such a system, interactions 
can be tuned in a broader range and would admit si- 
multaneous change of several pseudopotentials that may 
be required for the realization of other non-Abelian and 
multi-component FQH states [sil. Finally, we note that 
the experiments on the tunneling density of states [slj 
can directly measure the values of the tuned pseudopo- 
tentials [33|. 
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metrized version, |MR). The bare Coulomb state and the 
dielectric state behave very similarly under the change of 
^L, the latter being somewhat more robust. As noted in 
Ref. the overlaps with |MR) on the torus are large 
(97% and better), which is again the case for our dielec- 
tric state as well. 

Another independent insight into the nature of the di- 
electric state is the entanglement spectrum [sOl- The 
multiplicities of the low-lying levels of the entanglement 
spectrum contain topological information about a FQH 
state; this information is protected by an entanglement 
gap [2J whose magnitude measures the topological sta- 
bility of a state. We calculate and compare the entan- 
glement spectra on the sphere for the i/ = 1/2 Coulomb 
n = 1 LL state and the dielectric eigenstate tuned to a 
region of high overlap with the MR Pfafhan (Fig. EKb)). 
Not only do the both spectra display the same counting 
as the MR state (up to some value of the subsystem A's 
projection of the angular momentum Lf, set by the finite 
size of the sphere), but the entanglement gap (indicated 
by arrows in Fig. [SJb)) is also somewhat enhanced with 
respect to the unscreened Coulomb interaction. 

Conclusions. In summary, we proposed a method to 
tune the electron interactions in graphene and its bi- 
layer. In this approach, interaction pseudopotentials can 
be varied in a broad interval and FQH gaps can be en- 
hanced several times or even reduced to zero, allowing for 
a more complete exploration of compressible and incom- 
pressible phases than can be attained in GaAs materials. 
The proposed method is expected to be very efficient in 
optimizing the Abelian FQH states that belong to the 
hierarchy series, where the variation of the charge gap 
is followed by a small change in the overlap that gen- 
erally remains very close to unity. Non-Abelian states, 
although expected to be stabilized, may require a more 
subtle approach with several dielectric plates of different 
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